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Abstract 

We present a theory for the electron-temperature dependence Tg^ of opti- 
cal second harmonic generation (SHG). Such an analysis is required to study 
the dynamics of metallic systems with many hot electrons not at equilib- 
rium with the lattice. Using a tight-binding theory for the nonlinear sus- 
ceptibility x*'^H^;^e/) ^-ncl the Presnel coefficients we present results for the 
SHG intensity I^'^\io,Tf,i) for a Cu monolayer. In the case of linear opti- 
cal response we find that the intensity will decrease monotonously for in- 
creasing Tej. In agreement with experiment we find a frequency range where 
,T(.i) may be enhanced or reduced depending on electron temperature. 
Note, x^^H'^i ^ei) rather than the Fresnel coefficients determines essentially 
the temperature dependence. Our theory yields also that SHG probes effects 
due to hot electrons more sensitively than linear optics. We also discuss the 
Tej-dependence of SHG for Au and Ag. 
73.20.At, 78.47.+P, 78.66.-w 
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I. INTRODUCTION 



Due to its surface sensitivity, the nonlinear optical response has become a powerful probe 
for investigating the electronic structure of surfaces, interfaces, thin films, and multilayers. 
Recently, a combination of linear and nonlinear experiments has been performed exploit- 
ing the time dependence of SHG. Thus, effects of hot electrons not at equilibrium with the 
lattice and their changes in time can be analyzed. This opens a new route to investigate the 
dynamics of the system during relaxation to the equilibrium state. Note, different electronic 
temperatures of hot electrons not at equilibrium with the lattice but among themselves re- 
sult from varying the light irradiation. If intense and short laser pulses in the range of 100 fs 
to 1 ps are used for the SHG experiment, only the electrons will quickly thermalize (even in 
the approximation of Fermi liquid theory), since the slow electron-phonon coupling does not 
come into play Only later, at times of the order of several picoseconds the electrons 

will heat up the lattice. Thus, at short times the temperature dependence of SHG is essen- 
tially due to the varying non-equilibrium electronic temperature, which may be considerably 
different from the equilibrium temperature at later times, when electrons and lattice are at 
equilibrium. 

In this paper, we present a theory for the time-resolved optical SHG response, in partic- 
ular its dependence on non-equilibrium electronic temperature. We present results for the 
electron temperature dependence of the SHG yield at noble metal surfaces using an elec- 
tronic theory for the nonlinear response. The metallic surface is modeled by a freestanding 
Cu monolayer. However, the model allows already to identify characteristic features of this 
time-dependent nonlinear optical response to be expected also for the surface of bulk Cu and 
for other metals like Ag and Au and to explain recent experimental results for Cu and Au . 
Our analysis allows to identify the essential origin of the electron-temperature dependence 
of the nonlinear optical response. Furthermore our studies explain the different behavior 
of linear and nonlinear optics on non-equilibrium electronic temperatures. In general, our 
theory is of interest for the dynamics of nonequilibrium electronic systems. 

In section II, we describe details of our theory. In section III, we present results of 
our calculation of the nonlinear response for several values of the electron temperature 
and compare with experiments. The conclusions in section IV point out possible general 
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features of our model calculation and explain why the important temperature dependence 
results from x^'^H^^T) and not from the Fresnel coefficients. In the appendix some details 
of our calculation are presented. 



II. THEORY 



We calculate the SHG intensity I^^\uj,T) using an electronic theory. Then the SHG 
yield for p-polarization within the electric dipole approximation is ^ 



I{p - SH) = I (2.-) n Eo(^) n MiFcX^LVJs + N^F,{x'^Ul + cos^ ^ 



c 



+ N'Fsxfltl sin' cp) |2 . (1) 
Here, (p denotes the angle of polarization of the incident light, fc,s the Fresnel coefficients and 



ts^p the linear transmission coefficients, n = Je{uj) and = J ei2uj) the complex indices 



of refraction at frequencies u and 2uj. Furthermore, fs = smO /n and /c = \/l — fs for the 



fundamental frequency uj, Fg = sin9/A^ and = \Jl — F^ for the doubled frequency, 
where 6 and G denote the angle of incidence and the angle of reflection of the SHG light, 
respectively. The linear transmission coefficients are given by 



2cos^^ 2cos6' 2 cos 6 , . 

" ncosO + fc ' ^ cose + nfc ' ^ " A^cosB + F, ' ^ ^ 

The amplitude Ap in Eq. (1) is 

" cose • 

Note, the nonlinear susceptibility tensor x^"^^ is material specific and so are the linear dielec- 
tric function e{uj) and the indices of refraction n and via the electronic bandstructure. 

The contributions to I^'^\uj,T) due to the Fresnel factors Fc, Fs, fc, and fs and the 
transmission coefficients Tp, tp, and ts and x^'^K^^T) are all temperature dependent. In 
both cases the temperature dependence arises from the Fermi functions f{E,T) which, 
due to the many hot electrons resulting from the light irradiation, have to be taken at 
considerably elevated electron temperatures. Note, we limit ourselves to the time regime, 
where the electrons have already thermalized, but electron-phonon coupling has not become 
really effective. By taking both temperature dependences into account it becomes possible to 



decide theoretically whether x*^^)(u;, T) or the Fresnel and transmission coefficients cause the 
essential temperature dependence of the SHG intensities. Since it is known from theory and 
also experiment that the x^zL tensor element dominates over x^xlx x^zxx 05 we restrict 
our calculation to this single element of the nonlinear susceptibility for the SHG yield. Then 
approximately 



I{p-SH) = I (2^-) n Eoico) n A,N'F,xflfX^os'^ \' . (3) 



For the calculation of e{uj,T) and xEU"-^;^) ? we use eigenvalues -Ei^/.o- from a simple 
quadratic, freestanding Cu monolayer. The nearest-neighbor-distance equals 3.61/v^ A, ac- 
cording to the nearest-neighbor-distance of Cu bulk. The bandstructure involves five (i-bands 
and four plane waves and has been obtained [Q within the combined interpolation scheme 
(CIS). The (i-bands are parameterized in terms of Fletcher- Wohlfahrt parameters P,p!0 



their values are given in Table 1. The parameters of the ci-bands were evaluated by a fit to 



an ab initio LAPW band-structure calculation by Krakauer, Posternak, and Freeman [11 



Furthermore, to get also the correct onset of interband transitions we shifted correspondingly 



the d states |jT2|. The energy eigenvalues were calculated for 1861 k-points in the irreducible 
part of the Brillouin zone (1/8 of the whole Brillouin zone). The resulting band structure 
shown in Fig. 5 should be thus representative for the surface layer of bulk Cu. 

The calculation of the dielectric function e{uj, T) is performed by including both the intra- 



and interband electronic transitions. Neglecting the g-dependence of e, we take [|T3|,|T4| 



uj{uj + i/ Tpi) ^^^^ Eki>^ - Eua ~huj + thai 

where M are the dipole matrix elements, f{E,T) are the Fermi functions, u is the funda- 
mental frequency, Q is the unit area, ai = 0.1 eV is the Lorentzian broadening, Tpi the 
Drude relaxation time and Upi = ATm^e^ /m* ^ is the plasma frequency with n^. the electron 
density and m* the effective electron mass. Both ujpi and Tpi are fitted to a dielectric func- 



tion using literature values [0 at low fundamental energies, where no interband transitions 
occur. We use fiujpi = 9 eV and Tpi/h = 8.5 eV^"*^. 

To simplify our calculation, we assume constant matrix elements M, which fit to the 
dielectric function e{uj,T = OK) WM. This approximation is reasonable because the It- 



dependence of the matrix elements is expected to become less important in two dimensions 



due to the shrinking of the d-band width for the reduced coordination number and also due 
to the occurrence of additional allowed optical transitions. Additionally, as will be shown 
later, the main contribution of the SHG intensity has its origin in d ^ d ^ s transitions, 
so the matrix elements give a simple prefactor to the sum and cancel, when only intensity 
differences are considered. Once e{uj, T) is calculated, it is straightforward to get the Fresnel 
factors and transmission factors from Eq. (2). 

For evaluating the second-order susceptibility, we employ the microscopic theory developed 
in Ref. |T^. Thus, the tensor element xEz is given by 
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Xzzzi^,T) = — il4 <^ \ . (5) 

Note, for experiments not detecting absolute intensities (as is usually the case), the prefactor 
in this formula is of no further importance for comparison with experimental data. As in the 
calculation of e{uj) , we neglect the k-dependence of the matrix elements, taking the values 
for Mz from the e(co')-fit. We only take into account the dominant xElj-tensor element and 
neglect the contributions of the other tensor elements, as pointed out earlier. This completes 
then the theory for the determination of the SHG intensity I^'^\uj, T) with the two essential 
inputs e{ijj^T) and x*^^^(co',T) . 



III. RESULTS 

The most important results of our calculations are presented in Figs. 1-4. For these 
results we used an electronic structure which is shown in Fig. 5 and which was obtained 
using inputs discussed in section II. In Fig. 1 we present results for the change 

^ ' J(2)(^,300ir) 
of the SHG yield as a function of frequency for different non-equilibrium electronic tempera- 
tures Tel- In particular, in the inset of Fig. 1 we show for the frequency range from 0.4 to 1.5 
eV, that J*^^-* {uj) may be reduced or enhanced due to increasing non-equilibrium temperature 
for the electrons. In order to demonstrate what causes essentially the dependence on the 
electronic temperature, namely x^'^\'^iT) or the Fresnel coefficients, we present in Fig. 2 
results for 
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and in Fig. 3 results for the difference 

|x(2)|2 J(2)(cj,300ir)/|x(2)(u;,300i^)|2 

Note, in Fig. 3 we present results for I^'^\uj,T) using a frequency and temperature inde- 
pendent susceptibility x^^^- The comparison of the results in Figs. 2 and 3 show clearly 
that I^'^\uj,T) results essentially from x*^^H'^>^) ■ Furthermore, from comparing the results 
in Figs. 1 and 3 together with those in 2 and 3, we note that the crossover behavior with 
respect to the temperature dependence is only present in the case of the nonlinear response 
I^'^\uj,T) , but not in the Fresnel and transmission coefficients. 

In Fig. 4 we present results for I^'^\ijj,T — 2000ir) as a function of frequency. The 
different curves (b), (c), (d), and (e) refer to cases where we artificially excluded certain 
optical transitions from contributing to x'^^H^j ^) • Thus, we are able to demonstrate which 
transitions are of particular significance for the temperature dependence of I^'^\uj, T) . While 
curve (a) still includes all transitions, in the case of curve (b) we excluded contributions to 
X^'^\oj,T) , where the initial electronic state is not in the d-band. In the case of curve (c) 
we excluded contributions to x^'^\'^^T) , where the initial state is not in the d-band or the 
intermediate state is in the c^-band. Curve (d) results when the initial and the intermediate 
states arc in the rf-band, and in curve (c) we excluded transitions, where the initial state is 
in the d-band or the intermediate state is not in the d-band. Similar results are obtained 
for 300 K and 6000 K. 

For a detailed analysis of the results presented in Figs. 1 - 4 we like to make the 
following remarks. Investigating the relative changes of SHG intensity with temperature 
and its contributions due to the Fresnel factors and the nonlinear susceptibility x^'^\^^T) , 
we observe in Figs. 1, 2, and 3, that the shape of the dependence on the fundamental 
frequency is given by the second-order susceptibility x^'^\'^ iT) ioY fundamental energies 
below 3 eV. In the range between 1.2 eV and 3 eV, the SHG intensity decreases with 
increasing temperature. Most important is the result that indeed an energy window at a 
fundamental energy of 1.1 eV exists, where a small temperature increase (300 K to 2000 
K) results in a SHG intensity increase, but a stronger temperature increase (300 K to 
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6000 K) results in a rapid SHG intensity decrease. From Figs. 2 and 3 we conclude that 
X^'^\uj,T) is responsible for this. This clarifies then the physical origin of the temperature 
dependence of I^'^\uj,T) observed on Cu polycrystalline surfaces |Q, where in 2 eV photon 
energy pump-probe SHG experiments on Cu surfaces for weak pump pulses an increase of 
the SHG intensity has been observed, but a decrease when the pump pulse becomes stronger 
and thus the electron temperature higher. 

Note, it is remarkable that our model bandstructure yields already such a fair agreement 
with the experimental results. This is so, since the temperature dependence of T) and 

/*^^^(co',T) is strongly governed by the correct position of the upper d band edge and also 
because the important features of the bandstructure are already simulated correctly by our 
model. Furthermore, SHG probes the surface layer only, so our bandstructure is a fair 
approximation for the experimental situation. In view of this, we may also use with proper 
changes the model in Fig. 5 to describe the SHG response of Au and Ag. For this we change 
the onset of rf-band transitions A to A = 3.95 eV in the case of Ag and to A = 2.3 eV in 
the case of Au, since this is expected to cause the most important differences, s. Fig. 5. 
The other energy levels shift accordingly. 

Since the "mismatch" between {E^ - Ep) and the photon energy huj in Au is by 230 meV 
larger than in Cu (Cu bulk: Ed - Ep = 2.15 eV), one may simulate in our calculation the 
Au surface by using approximately the SHG spectrum of the Cu monolayer, but considering 
a photon energy which is 230 meV below the respective energy value for the Cu monolayer. 
While for Cu we found the crossover behavior at 1.1 eV, we have to take for Au a fundamental 
photon energy of 0.87 eV. Then our results presented in Fig. 1 yield an increase of SHG 
for all temperatures. The increase for 6000 K is higher than that for 2000 K. These results 
seem in reasonable agreement with experiments for a Au polycrystalline surface, where a 
monotonic increase of the SHG yield for both temperatures has been observed Q]. Thus, 
our electronic theory is able to explain the ultrafast electronic relaxation process on both 
Cu and Au surfaces. From our direct calculation of the dielectric functions (see Fig. 6) we 
find that the temperature dependence of J*^^-* is mainly caused by x^'^K^^ T) . At a frequency 
hu = 0.9 eV we find that the x^'^K^^T) -dependence on the electron temperature begins 



to saturate for higher temperature, as is observed in the experiment |^,p!8|. Note, the SHG 
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response of Ag may be modeled similarly as has been described for Au. Further remarks on 
this are given later. 

Regarding the origin of the temperature dependence of the SHG contribution (on transi- 
tions), it is necessary to investigate the contribution of the rf-band electrons to ^^}jyu)^ T) and 
its influence on the SHG yield in some detail. Results of this analysis are shown in Fig. 4. 
To interpret the results, one has to be careful not to neglect the interferences between the 
various terms in Eq. (5). Therefore, calculating xfJzi^, T) and the SHG intensity for just a 
few transitions, neglecting all others, is of little value, since then the interferences between 
the various complex quantities are almost completely neglected. Thus, we calculated the 
SHG intensity with a xf)z{'^i ^) ' Authors: T. A. Luce, W. H"ubner, and K. H. Bennemann 

where just some transitions with a specific combination of the three states necessary 
for a nonlinear transition were neglected. States were identified as d-states if their energy 
eigenvalue was between -5.25 and -1.45 eV below Ep. Note, here we neglected the symmetry 
of the states. Apparently, we find that the initial state is almost always a d state, since both 
intensities nearly coincide. The SHG intensity /*^^)(a;, T) from transitions d —>■ {s,p) — > 
{d,s,p), where the initial state is in the d-band and the intermediate state is not (curve 
(c)), nearly vanishes. The ratio of the maximum values of the SHG intensity and the total 
SHG intensity I^^\d {s,p) {d, s,p))/ I^'^\total) is equal to 1/90. Of course, both 
intensities are taken at the same frequency. This shows that most transitions contributing 
to the SHG intensity have the initial and intermediate state in the d-band. This is confirmed 
by the respective dotted curve (d) . Due to interference of the various interband transitions 
in xEU*^)^) ) this curve and the total SHG intensity yield differ at the maximum by a 
factor of 3. Comparing this ratio I^'^\d ^ d ^ {d, s,p)) / I^'^\total) = 3 with the ratio 
I^'^\d {s,p) — ^ (d, s,p))/ I^^^total), it is reasonable to assume that only the transitions 
with the initial and the intermediate state in the d-band contribute significantly to the SHG 
yield. 

From curve (e) in Fig. 4 it can be seen that considerable SHG intensity is created by 

(s.p) ^ d {s,p,d) transitions. This is readily understood if one notes that all terms in 

the sum of x^^Jzi'^^ ^1- consist of differences 

1 f f{E,. T) f{E2. T) fjE,. T) f{E,, T) \ 

E3- El- 2T1UJ + 2ihai [Ea- E2-hu; + ihai E2 - Ei - ?iuj + ihai j ' ^ ^ 
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Here, Ei, E2 and E3 denote the initial, intermediate, and final state in the considered 
transition, respectively. This difference will give a contribution even if one of the terms in 
curly brackets is zero, so the intensity I^'^\s,p d ^ s,p, d) mainly results from the first 
term. 

Compared with these transitions, all other combinations of states contribute only a 
small amount to the SHG. Even at low energies {fiw ~ 1 eV), the probabihty that all 
three states of a nonlinear transition are d-states is rather small. Since the main con- 
tribution to the SHG intensity results from transitions with the two lower states in the 
(i-band, nearly all contributing terms have the corresponding dipole matrix element product 
{d\z\d){d\z\{s,p)){{s,p)\z\d). This can be handled as a prefactor. In view of this result, our 
simplification of using one value for all dipole matrix elements is justified. 

In order to avoid confusion, one should note that our calculation of 
I^'^\uj,T) demonstrates that the mismatch energy S = Ep — E^ — hcu is crucial for the 
T-dependence of the SHG yield and gives for a fixed photon energy hu the following results: 
First, if (5 < Scu, where Scu refers to the value where the refiectivity of Cu saturates, the 
temperature dependence of the Presnel factors and transmission coefficients is sufficient to 
describe I^'^\uj,T) . In the range 5cu < S < 5au, I^'^\uj,T) is essentially due to x^^K^,T) , 
and for 6 > 6au, I^'^\^,T) results from x*^^)(c<j,T) for both Cu and An. This is consistent 
with experimental observations regarding the dependence of the refiectivity on 5. 

In summary, our results for I^'^\u,T) suggest that the temperature dependence is es- 
sentially due to x*^^)(a;,T') . In particular, x(^)(a;,T) is responsible for the nonmonotonous 
temperature dependence I^'^\uj,T) . 

IV. CONCLUSIONS 

We calculated the dependence of the SHG yield on electron temperatures for electrons not 
at equilibrium with the lattice, caused by light irradiation for a Cu surface and for a Ag and 
a Au surface. We find that for energies below 3.2 eV the temperature effects result mainly 
from x^^) (a;, T) and not from the Presnel coefficients. In particular, our calculation yield that 
in the frequency range between 0.8 and 1.1 eV, the SHG intensity even may increase with 
increasing electron temperature. This is the case for Cu at low light intensities and for Au for 
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all light intensities that cause no damage. Note, for Au we use the same bandstructure and 
smaller light frequencies in order to study the same electronic transitions as for Cu. Clearly, 
/*^^^(co',T) depends on the position of the cZ-band with respect to the Fermi-energy. The 
interesting effect that mainly x'^^)(ti;,T) causes the temperature dependence comes about 
due to extra two photon absorption processes making use of the high (i-band density of 
states and occurring only at elevated electronic temperatures. Note, however, that only due 
to the interference of these transitions with others a considerable temperature dependence 
comes about. The situation is illustrated in Fig. 7. Note, without the effects due to the 
indicated transitions, the nonmonotonic temperature dependence would be absent. The 
figure illustrates also the difference between the temperature dependence in nonlinear and 
linear response, where for the latter only (i-states as initial states are possible. 

As already remarked the essential features of our results should be valid independent 
of our model calculations. This is supported by the fact that we already find such a fair 
agreement with the experiments by Hohlfeld, Conrad and Matthias Q]. Our theory shows 
clearly how I^'^\uj,T) depends on the electronic structure and thus what can be expected 
for other noble metals and transition metals. For example, for Ag {Ep — Eii = 3.98 eV 
for bulk Ag) we would expect, using the previously discussed argument, a similar behavior 
for the linear and nonlinear response in the frequency range from 0.8 to 1.1 eV in view of 
Fig. 7, while in the frequency range 4 eV a decrease of the SHG yield due to the very 
small negative mismatch is expected. For Ag and for this mismatch the SHG intensity 
I^'^\uj,T) should always decrease as a function of T. Also for Fe, in addition to the changes 
due to the electronic structure (position of the d-band) we expect interesting non-equilibrium 
temperature effects due to the additional dependence of the magnetization M on T^i, since 
J*^^^ depends on T^i and M(Te/). 

Regarding the time dependence of our results, for larger times the electron-phonon cou- 
pling will become more effective and cause the decrease of the electronic temperature via 
an energy transfer from the hot electrons to the lattice. The resulting equilibrium temper- 
ature will be only a few hundred degrees higher than before the light irradiation. Then, 
/(2)(cu, T) will be similar to I^^\uj, 300K), which is given in the appendix. Such temperature 
dependences should be compared with those obtained for systems due to usual equilibrium 
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thermodynamical effects. It is straightforward to extend our calculation to thicker films 
consisting of several atomic layers. Such calculations are in progress []TU| . 



As a resume, our studies show that SHG can be used to study the dynamics of excited 
hot electrons in crystals. Note, in this paper we considered the time window where the 
electrons are far from equilibrium and have their own temperature different from the lattice 
temperature. Obviously, using SHG for studying the time evolution of the electronic system 
far from equilibrium offers new perspectives, in particular for studying magnetism. 
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APPENDIX 

In this appendix we present some details of our calculations. In Fig. 5 we show results 
for our band-structure calculation, evaluated with the parameters fitted to the ab initio 
bandstructure. Our CIS calculation is in good agreement with the ab initio calculations. 
The resulting density of states (DOS) in the inset of Fig. 6 gives a (i-band edge approximately 
1.45 eV below the Fermi-energy Ef. The total c/-band width is 3.05 eV. (i-band transitions 
with lowest energy are possible at energies around 2.1 eV (near the X-point). The constant 
values of the matrix elements are = 5.4 * 10^^° m. The numerical value was obtained by 



fitting the maximum values of the dielectric function to literature values [|T^. From atomic 
orbitals, one would expect values about 0.01 * 10^^° m for the {p\z\d) matrix elements [jl6|. 
This shows that the atomic orbitals give the correct symmetry of the wave functions, but the 
absolute values differ considerably due to the periodicity of the monolayer wave functions 
which delocalizes the d-electrons. 

The real and imaginary parts of the dielectric function for the Cu monolayer as a function 
of the incident photon energy are plotted in Fig. 6 for three different electronic temperatures 
(300 K, 2000 K, and 6000 K). These temperatures are chosen as realistic estimates of the 
electron temperature caused by experimental pump-probe pulses. 

The influence of the plasma frequency ujpi at low frequencies leads to a strong drop 
(increase) of the real (imaginary) part of e{uj) upon lowering the frequency uj. A variation of 
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the Drude relaxation time Tpi causes only minor changes • If ujpi is reduced, the minimum 
of Im {e} increases in value and moves to lower energies. The increase of Im {e} at 2 eV 
is caused by transitions from the d-hand to unoccupied states at Ef at this energy. Due 
to the Lorentzian broadening introduced by the calculation, the increase of Im {e} is not 
as steep as expected from optical measurements. Rising the electron temperature results 
in a flattening of the curves for Re {e} and Im {e}. The minimum of Im {e} becomes less 
pronounced for higher temperatures. Since there are more vacant states below Ef a.t higher 
temperatures, interband transitions occur at lower energies, thus weakening the rise of Im 
{e} at 2 eV. It is to be expected that the main temperature effects are observable at photon 
energies u with either hu ^ Ef — E^ or 2huj ^ Ef — E^, depending on the influence of the 
Fresnel factors or the nonlinear susceptibility on the temperature dependence. Then the 
transitions are strongly influenced by the effects which a temperature increase induces on 
the occupation of states at Ef. 

In order to trace back the origin of this effect we separately calculate the contribution of 
the Fresnel factors and nonlinear contributions to the SHG yield. To obtain the contributions 
of the Fresnel factors we calculate the SHG intensity I{2uj) from Eq. (3) with x^^l set to 
unity, thus including the interaction of the incident light and the frequency-doubled light 
with the monolayer. For the nonlinear contributions we compute x^^zzi^^'^) ■ Inserting the 
dielectric function from Eq. (4) and xizlli^? ^) front Eq. (5) in Eq. (3), in Fig. 8 we 
show the calculated SHG intensity I^'^\uj,T) for the three different electronic temperatures 
(300 K, 2000 K, 6000 K). The SHG intensity rises strongly at photon energies around 2 
eV and peaks at approximately 2.8 eV. Then it drops without exhibiting further peaks. 
For increasing temperature, the intensity maximum decreases and moves to slightly higher 
photon energies. From Fig. 9 we find a similar shape of the xEll^-^; T) -plots, with maxima 
shifted to photon energies around 1.9 eV. From Eq. (5) it becomes clear that even at low 
energies in the range of 1 eV contributions to xEU'^;^) ^-^e to be expected due to 2huj- 
resonances of the form 

1 

-£^k+2q,«"<7 - -£^k/a - 2^0; + i2hai 
The maximum decreases with rising temperature, but its position remains at 1.9 eV. Fig. 10 
displays the contributions of the Fresnel- and transmission factors to the SHG intensity. We 
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find a strong increase at energies larger than 2 eV, due to the interband transition threshold 
in this energy range. Decreasing the plasma frequency increases the SHG intensity and shifts 
the maximum to lower frequencies (e.g. for ojpi = 2.2 eV we find a maximum at 1.1 eV, 
the intensity amplitude there is about 1000 times larger). The plasma frequency influences 
both Fresnel factors f{uj) and F{2uj). However since it is necessary that both F{2uj) and 
f{uj) are nonvanishing to make SHG possible, in the Presnel factor contributions of the SHG 
yield mainly the effect of f{uj) is visible, because F{2u;) ^ for photon energies greater than 
energies where is nonvanishing. 

Analyzing the contributions of the linear and the nonlinear response to the SHG intensity, 
we find that the overall shape of the intensity is dominated by the xizU*^, 7") - contribution, 
if compared with the corresponding ■^^}JyLij^T') -plots in Fig. 9. The apparent shift of the 
xi;!^'^' ^) -plot to higher photon energy results from the strong increase of the Presnel factor 
contribution to the SHG in this energy range, as can be seen from Fig. 10. 
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TABLES 

TABLE I. Fletcher - Wohlfahrt parameters fitted to ab initio calcu- 
lations for the Cu monolayer. (E'er On-site energy, Af Overlap integrals 
of the d-orbitals) 



parameter 


value [eV] 


Eq 


-2.904 


Ai 


0.106 


A2 


0.088 


A3 


0.147 


Aa 


0.220 


As 


0.170 


Ae 


0.200 
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FIGURES 

FIG. 1. Change of the SHG intensity A/(2)(a;) = ^^''^'^yTa'j'l^sooK)^""'^^ *° * of electron 
temperature as a function of the incident photon energy hto in percent for T^i = 2000i^ and 
Tel = 6000iir. The inset shows the energy range between 0.4 and 1.5 eV at an enlarged abscissa 
scale. 

FIG. 2. Change of the second order susceptibility Ax^'^\u;,T) = ^^^^^'^''^^i^J^^QK)^^^^ 

(2) 

a rise of electron temperature as a function of the incident photon energy Tiu in percent. Xzzz is 
very similar to the SHG yield for energies below 3.2 eV, indicating that the temperature effects 
for these fundamental energies on the SHG intensity are caused by xizz- The inset displays this 
similarity for incident energies between 0.4 and 1.5 eV. 

FIG. 3. Results for the relative change of the SHG yield 

contrast to the results shown in Fig. 1, the frequency and temperature dependence of x^^H'^;^) is 
neglected. The inset shows at an enlarged abscissa scale the energy range between 0.4 and 1.5 eV. 
These results demonstrate the importance of the temperature dependence of x^'^\^iT) ■ 

FIG. 4. Contributions of the different electronic transition combinations to the SHG inten- 
sity as function of the incident photon energy hu (arbitrary units). Note, the contributions due 
to transitions with the incident state a d-band (curve (b)) nearly coincides with the total SHG 
intensity. 

FIG. 5. Calculated bandstructure of a Cu monolayer using the Combined Interpolation Scheme. 
The d band edge is 1.45 eV below Ep- A is the interband transition onset energy. The symmetry 
of the bands is indicated by the s and d labels. 

FIG. 6. Real and imaginary parts of the dielectric function as a function of the incident photon 
energy Tiuj for three nonequilibrium electronic temperatures (300 K, 2000 K, 6000 K), calculated 
from Eq. (4). The inset shows the DOS of the Cu monolayer bandstructure. 



16 



FIG. 7. Illustration of important electronic transitions generating SH from noble metals. Note, 
both initial and intermediate state belong to the d-band. Actually, by neglecting such transitions 
in our calculation, we find a drastically reduced I^'^\iv,T) . Since in the case of the linear optical 

response, at best only the initial state can be a d state, noncquilibrium temperature effects are 
drastically reduced and less interesting with regards to probing the electronic structure. 

FIG. 8. Calculated total SHG intensity of a Cu monolayer for three different electronic temper- 
atures (300 K, 2000 K, 6000 K) as a function of the incident photon energy hu) (arbitrary units). 
This quantity has been used for the calculation of the intensity changes. 

FIG. 9. Calculated total susceptibility Xzzz for three different temperatures (300 K, 2000 K, 
6000 K) as a function of the incident photon energy huj (arbitrary units) . These results have been 
used for the calculation of the susceptibility changes given in Fig. 2 . 

FIG. 10. Temperature dependence of the SHG intensity resulting from the temperature depen- 

(2) 

dence of the Fresnel factors and putting Xzzz{i^,T) equal to unity in Eq. (3). These results have 
been used to calculate the differences in Fig. 3. 
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